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Abstract 

CsJ ' In the framework of perturbative QCD and a model for the production of charmed 

I hadrons the structure function F| is calculated and compared with the experimental data 

I of HI and ZEUS Collaborations. We show that the spectator mechanism of the D*-m.eson 

pi ^1 production independent of a hadronic remnant is valid at pT > 10 GeV, only. We find 

I that the evolution of Q'^) versus the virtuality can be neglected in the kinematic 

' region of HERA. 
(N ■ 

> '■ 1 Introduction 
00 

\o _ 

^ i As well known, the DIS data in a broad kinematic region can be precisely described by a 

O . set of universal partonic distributions obeyed equations of the DGLAP-evolution at large 

Q ■ virtualities Q^. Some difficulties appear when one tries to consider the heavy quark contribution 

^ ' into the function F2 in the DGALP-technique. The problem is caused by an accurate account 



X 



of the kinematic region Q ~ 4m^, where, on one hand, we deal with high virtualities providing 
Q I the applicability of pQCD, and, on the other hand, the mass effects should be appropriately 
treated and the heavy quarks cannot be considered as massless. Some attempts to take into 
account the heavy quarks contribution near the threshold and to match this contribution with 
the general DGLAP equations at large ^ Am'^ are presented in ^. In this work we 
consider c-quark production by making use of the notations in , where the attempt to match 
two different approaches for large and ~ 4m^ was done. 

The c-quark contribution into the structure function in the deep inelastic scattering can be 
obtained by convoluting the partonic distribution with the coefficient functions: 

^F^{x,Q') = Cg{Q\fi')®fgifi')[x] + 
+C,(g^/i2)®/e(/x^)[x] + 

(1) 

where denote the partonic densities in the proton, the symbol ® corresponds to the convo- 
lution over the variable x: 

dz 



b[x] = I —a{z)h (-) 



At small iJ? < fil the charm contents is /c(x,/i^) = 0, and the structure function is 
completely described by the photon-gluon fusion Y9 ~^ cc, that gives Cg{Q'^,iJ?) ® fgifJ^'^)- 
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The second term in (|I|) corresponds to the process of virtual photon scattering on the c-quark 
from the sea of the initial hadron. At large the process 7*c eg becomes dominant. The 
process 7*g — > ccq is usually not taken into account because its contribution is suppressed by 
the additional power of a^, however below we show that this process can play an essential role 
in the region of moderate pt and Q^. 

In addition to the problem of matching the different approaches in the kinematic regions 
of ~ ml and 3> 4m^, another important problem is connected to the comparison of 
theoretical predictions with the experimental data. So, F"^ is reconstructed in the ZEUS and HI 
experiments over the D*^-meson production data obtained in two decay modes: D* — > K27r, 
D* KAtt. The D*-meson production cross section is measured in the restricted kinematic 
region: at ZEUS one puts 1 < < 600 GeV^, \r]{D*)\ < 1.5, 1.5 < p{D*) < 15 GeV, 0.02 < 
y < 0.7, while the HI Collaboration uses the similar kinematic region: 1 < < 100 GeV, 
\t]{D*)\ < 1.5, 1.5 < Pt{.D*) < 10 GeV, 0.05 < y < 0.7. The cross section in the whole 
kinematic region is reconstructed under some model assumptions of charm production. One 
can conclude that the really measured quantities are the spectra of the charmed mesons in the 
restricted kinematic region, but the c-quark ones. That is why one has to define a procedure of 
the charmed quark hadronization in order to compare the calculated function with 
the experiment. 

A procedure commonly used is based on the factorization theorem applied for the momentum 
spectrum of D*-meson: 

- — = — ®D{z), (2) 
dpT dkx 

where D{z) is the fragmentation function of c — » D*, and z = Px/kT- The form of D{z) is given 
by some kind of phenomenological anzatze [Q. Equation (Q) should be a good approximation 
at large pt and kx-, where the factorization theorem is valid. This theorem cannot be applied 
to the whole kinematic region since this could generate essential errors in the relation between 
the D*-m.eson and c-quark spectra. The drawback of is clearly seen in the HI and ZEUS 
data: there is a visible excess of events in the region of pseudorapidity distribution towards 
the fragmentation of proton-beam. One can see the same features in the distribution over 
the variable zj:,* = {Pp£,*)/{Pq), where P, po*-, Q are four-momenta of the proton, Z}*-meson 
and virtual photon, correspondingly. The discrepancy between next-to-leading-order (NLO) 
QCD prediction under the hadronization mechanism of @ and the experimental data is more 
essential at large rjjy* and zd*. 

Usually one tries to improve the data description by using some additional model assump- 
tions incorporated in Monte-Carlo codes of event generator 0]. 



2 The charm production model 

We will discuss the model where the attempt to take into account the hadronization from 
the very beginning is done. Let us consider the all pQCD 0(aa^)-diagrams, describing the 
perturbative production of the cg-pair with the quantum numbers of the appropriate charm 
mesons. 



2 



A complete set of diagrams is presented in Fig. 1. In this chapter our consideration is 
restricted by the case of D*-meson production in the photon-gluon interaction '-)*g ^ D* + 
c + q (see equation (|lD). One needs the additional parameter (O) to take into account the 
contribution by the simultaneous production of c - and g-quarks. This parameter describes the 
fusion of c-quark and its light co-mover q into (cg)-system, and that is proportional to the D*- 
meson wave function squared. This parameter can be extracted from the experimentally known 
value of the fragmentation probability W of c ^ D* at large transverse momenta, W = 0.23 
i- 

The (cg)-system can be in two color states: octet or singlet. In our previous analysis of the 
data on the charm photoproduction and electroproduction [0 we considered the production of 
color octet and color singlet states independently. So, we used two fusion parameters: 
and (0(8)), correspondingly. The ratio of these values is close to unit. At small pr the color- 
singlet contribution into the cross section as well as the color-octet one behave as and 
they are approximately equal to each other. At large p-r both contributions behave like (1/p^), 
and the singlet term is dominant. The color-singlet dominance takes place due to the color- 
coefficient ratio. The both contributions at large pt come to the fragmentation regime. In 
this regime the octet and singlet have a similar behavior, and the production process can be 
characterized by the only parameter: 

(0^ff) = (0(i)) + ^(0(8)) 

. A more detail description of the model is represented in the Appendix. 

One of essential features of this model consists of obeying the factorization theorem at large 
transverse momenta. In other words, at large pt the formula of (0) describes the process with 
a good approximation, and the process looks like a fragmentation. At small momenta the 
fragmentation form of the charm production is strongly broken by additional terms ~ 
contrast to the fragmentation term, which behaves like 1/px cis follows from (0). 

One can see from Fig. 2 how the factorization theorem works at large pt for both the 
singlet contribution and the octet one. The contributions of complete diagram-set (Fig. 1) in 
comparison with the contribution of the fragmentation diagram one (as in (^) are presented 
in the figure for both the singlet and octet. 

The calculations have been done at rather large value of s^g = {p^ + Pg)'^ = 200^GeV^ in 
order to reach the factorization regime (|]) . One can see from the figure that there is the region 
of Pt, where the fragmentation mechanism is not valid, and the contribution of perturbative 
recombination is dominant. The presence of the recombination contribution is due to the gauge 
invariance of QCD. The p^-dependence of this contribution has the additional factor l/p^, that 
is why the recombination vanishes at large pt, and the D*-meson production is described by 
the factorization formula of (^. 

3 The D*-meson electroproduction 

As was mentioned above, the model under consideration has been used to describe the charm 
photoproduction data of ZEUS Collaboration [Q. The fragmentation probability W{c D*) 
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was extracted from the experimental data on the D*-meson production in e^e~-annihilation 
1^. The value of factorization scale for the operators (0(i^8)) has been fixed at the Z)*-meson 
mass: fip = mo*- The light quark mass equals rriq = 0.3 GeV, and the c -quark mass equals 
m, = 1.5 GeV. 

The factorization probability W{c D*) = 0.23, while (O^^) = 0.25 GeV^. The ratio 
between the octet and singlet operators has been chosen equal to (0(8))/(0(i)) = 1.3. It is 
worth to mention that taking into account the color-octet contribution essentially improves 
the description of the experimental data, especially, in the description of the pseudorapidity 
distribution. In this distribution the octet enforces the production into the forward semisphere 
(toward the direction of the initial proton) and improves the data description in comparison 
with the NLO prediction 0]. 

The calculations of cross section for the D*-meson electroproduction have been made by 
practically the same manner as for the photoproduction. All we need is to replace the density- 
matrix of real photon by the density-matrix of virtual photon. So, after averaging over the 
photon and electron polarizations the D*-meson elecroproduction amplitude squared has the 
following from: 

PP _|_ p t* _ Sln^j 

where \A\'^ is the D*-meson electroproduction amplitude squared, Mj is the D*-meson pho- 
toproduction amplitude {i is the Lorentz index of the photon polarization), /ciand /c2 are the 
initial and final positron momenta, and = —{ki — A;2)^- The calculation results are presented 
in Figs. 3 and 4 in comparison with the experimental data of the HI and ZEUS Collaborations, 
correspondingly. The calculations have been done with the same cuts as in the appropriate 
experiments. The singlet contribution is plotted by the dashed curve, the octet contribution is 
plotted by the dotted curve; the solid curve corresponds to the sum of these contributions. 

We have used two values of the scale for the running coupling constant of QCD in our cal- 
culations: fiR = \/mjy, + (upper curve) and hr = a/ 4m|,* + Q"^ (lower curve) and the scale 
value /iF = A/ml,, -|- Q'^ for the CTEQ4 parameterization of the structure function. As one can 
see from the figures the model prediction is in agreement with the experimental data. Such the 
good description is achieved due to taking into account the octet contribution. This contribu- 
tion improves the distributions over the pseudorapidity f]{D*) (as in the photoproduction); over 
the variable z{D*) measured in the HI experiment, and over the variable x{D*) = \'p'{D*)\/W 
measured in the ZEUS experiment {p*{D*) is the D*-meson 3-momentum in the c.m.s. of initial 
virtual photon and proton). It is worth to mention another feature of the octet contribution. 
It is essential at small pt and becomes negligible at large ones, where it is suppressed by the 
color-factor 1/8. 

Thus, in the framework of the model under consideration we have achieved a good descrip- 
tion for the ZEUS charm-photoproduction data P] as well as for the charm-electroproduction 
data of HI and ZEUS Collaborations |Q. This circumstance allows us to suppose that the 
extrapolation of the experimental data into the total kinematic region with the help of our 
model is rather reliable. It means, that we can calculate the total cross section production and, 
therefore, extract 
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4 The structure function F2 

The charm contribution F2 into the structure function F2 is defined by the doubly differential 
cross-section of the charm-production as follows: 

{[l + d-l-n^l-i/W..)} (4) 

Generally one neglects the contribution of the longitudinal component because of its suppres- 
sion. As was mentioned in the Introduction, F| is reconstructed on the base of the experimental 
data for the D*-meson production in the experimentally available kinematic region. The ob- 
served production cross section is extrapolated into the total kinematic region in the framework 
of some model. Thus, it is clear that the F2 value depends on a model. In the framework of 
our model we have an opportunity to calculate the cross section in the total kinematic region, 
determine , and compare it with the results of ZEUS and HI Collaborations. 

The experimental dependence of -F| on x for different values is shown in Fig. 5 by dots 
(the HI Collaboration). These data were extracted by the extrapolation of the experimentally 
observed cross section into the total kinematic region on the base of NLO calculations and 
Monte-Carlo programs taking into account the hadronization. 

The curves in this figure correspond to our model predictions. -F| has been calculated 
according to formula The falls on the distribution tails appear because of the phase-space 
borders for the given value of the ep -interaction energy and chosen values for the quark masses. 

The ZEUS experimental data on F2 are shown in Fig. 6 in comparison with the predictions 
of our model . 



5 The perturbative recombination 

One can see that the operator expansion of ([l|) contains not only the term interpreted as photon- 
gluon production of charm, but also that of describing the photon-quark production of charm. 
Generally one neglects this term due to the additional factor of a^. 

The diagrams, which correspond to the photon-quark term in our approach, are shown in 
Fig. 7. For such kind of diagrams the only difference between the color-octet production cross 
section and the color-singlet one is due to the overall color factor of 1/8. That is why we can 
restrict ourselves by considering the singlet production, only. 

The recent analysis demonstrates that at large transverse momenta pt the photon-quark 
production contribution is suppressed by an additional factor l/p"^. At small pr the suppression 
is absent, and the differential cross section of the (eg) pair production at the angle = (toward 
the direction of initial quark q) has a large numerical coefficient in comparison with the c-quark 
production in the 5f7-interaction 0: 

da{'y + q ^ (cq) + c) 2567r 

da{^ + g^cc) - ^ ^ 

Thus, the smallness of as in the photon-quark production cross section is compensated 
by large numerical coefficient. To the same moment, the production at the angle B = vr is 



5 



suppressed by the additional power of energy. Another circumstance is essential, too: this 
contribution slightly depends on the light quark mass and does not vanish at — > 0. 

In paper the cg-pair production cross section in the interaction of the photon and valence 
quark from the initial proton has been calculated to evaluate the /D~ and D^/D^ asym- 
metries. The predicted value of asymmetry in the kinematic regions researched by the ZEUS 
and HI Collaborations is about 2-3%, that is in the same order of magnitude as experimental 
errors. At low energies of 7p-interactions, the role of the photon-quark production becomes 
essential and this contribution yields the asymmetry prediction which is in a good agreement 
with the experimental data 0. If the production asymmetry is due to the perturbative recom- 
bination indeed, then the asymmetry decreases with the pt increase, because the perturbative 
recombination has the additional l/pr factor in comparison with the leading contributions. On 
the other hand, if the asymmetry is due to the interaction between the c-quark and the valence 
quark from the initial hadron, then the interaction between the c-quark and the light quark 
from the initial hadron sea exists. Such the contribution has been calculated in the framework 
of our model, and it surprisingly looks like the octet contribution in the distribution shape as 
well as in the absolute value. It is worth to mention that the quark-photon contribution does 
not contain an additional normalization factor, which the octet one contains. In Fig. 7 both 
the octet and quark-photon contributions into the D*-meson production are presented for the 
kinematic region investigated by the HI Collaboration. One can see that the distributions over 
logio(a;) and px are practically the same for the whole investigated range. One can see the 
only small difference in the normalization. The Q^-distributions at < 10 GeV^ have prac- 
tically the qualitatively similar behavior, too. The differences between the distributions over 
W, 1]{D*) and z{D*) are more essential for the production mechanisms compared. However, 
these distributions have qualitatively analogous behavior. The singlet cg-pair yield in 'yg- and 
7g-interactions is presented in Fig. 8 for the kinematic conditions of HI experiment. It is clear 
from this figure that at not very large the sum of the 'yg- and 7g-contributions into the 
singlet describes the data so good as the sum of 75f-contributions into the singlet and octet. 

Therefore, the singlet cg-pair production mechanism is enough to describe the charm pho- 
toproduction and electroproduction data in the HERA experiments. The photon-gluon con- 
tribution as well as the quark-photon one play essential role for the singlet cg-pair production 
at the HERA interaction energy. The serious arguments to take into account the sea quark 
contribution into the charm production was presented in paper ||^, where the charm production 
asymmetry has been successfully described in the framework of perturbative recombination in 
the E687 and E691 experiments. 

6 Conclusions 

The model under consideration is based, at first, on the heavy-quark production in the pertur- 
bative theory, and second, it uses the nonperturbative model of quarks fusion into the hadron. 
This model allows us to describe the existing data on charm photoproduction and electropro- 
duction in the total kinematic region. 

Let us itemize the main features of the model predictions: 
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1. In the region of large transverse momenta of the cg-pair our model predictions coincide 



mentum spectrum of cg-pair is calculated by convolving the heavy quark spectrum with 
the fragmentation function. 

2. In the region of small transverse momenta and small the main contribution into the 
inclusive spectrum is due to the recombination diagrams, that depends on pt as in 
contrast to the fragmentation ones, which depend on pt as The contributions of 
(77- and g7-interactions into the charm production are comparable. 

3. It is not necessary to include the term Cc®c of F2 (the second term in expansion(P) into 
the data description at the HERA energies. According to the paper the term Cc® c 
becomes essential at large Q^. 

Therefore our main conclusion is formulated as follows: the spectator character of the Z}*-meson 
production, which is independent of the flavor of the initial hadron remnant, becomes dominant 
at > 10 GeV. At small transverse momenta the essential part of the cross section is due to 
the interaction with the initial hadron remnant. The interaction with the valence quark of the 
remnant explains the experimentally observed flavor asymmetry in the charm yield. 
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Appendix 

A. The hard production of four quarks in the photon-gluon subpro- 

cess 

As was mentioned above, in the model under consideration one supposes that the both valence 
quarks, the heavy and light ones, are produced in the hard process. 

In the framework of the tree level approach to the subprocess g{pg) + 7(^7) c(gc) + 
c(gc) + d{qq) + q{qq) 24 Feynman diagrams contribute (see Fig. 1; he bold line corresponds to 
the c-quark, the thin line corresponds to the (i-quark). In this paper the amplitude calculation 
has been done by straightforward multiplying of 7-matrices, spinors and polarization vectors. 

Let us introduce the following notations: 



with the predictions of the factorization model in form (0), or in other words, the mo- 



22 2 
(lc=<l-c= rric 




Pgj Pg ~l~ P-yi 



Qqq Qq Qq 



■g,cc — 



Pg + qc 



'g,qq — + Qqq- 



(6) 
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The three-gluon vertex has the form 

+ [{2pb + Pa)-ea)€l 

^'^{Pa,Pb,ea,eb) = T^{Pa,Pb,ea,eb)/{pa+ Pbf- (7) 



The bi-quark currents are defined as follows: 

•Jcc 2_ ' 99 „2_ ■ V / 

Also we need to define the auxiliary spinors 

-/ iQciq) ~ Pg(-y) + mc(g)) (^9(7) " ^5(9) + "^g(g) ) „ 

Uc{q)g{y) = U{Qc{q))eg(y)J- —3—, ^^5(5)9(7) = 7" " —^^3(^)^(5^), 



u{qc)Jqq{qc- qqq + rric) _ u{qq) Jcc{qq - Qcc + m. 



^"''^'^ {qc-qqqY-ml ' {qq - qcc)'^ - rn? 
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_ {qqq- qc + mc)JqqV{qc) _ {qcc - Qq + fUg) JcgV jqg) , . 

{qqq-q.Y-K ' fe,-g,)^-m^ ' 

where eg and are the polarization vectors of the gluon and photon, correspondingly. The 
matrix element squared is summed over the following ortonormalized states of gluon: 

6' = (0,1,0,0), e" = (0,0,1,0). 

It is worth to mention that e'^ = e"^ = — 1, and ]? ■ e = 0, where p is gluon momentum. 

In the case of deep inelastic production the photon is off mass shell, and one need to replace 
the matrix €2 in (|^) by the matrix 7*. Thus the matrix element will have a free Lorentz index 
for convolving with the photon density in accordance with formula (|^). 

We use the following index definitions: 

the upper indices jg designate the color state of gluon; 

the low index designates the color state of c-quark; 

the low index designates the color state of c-quark; 

the low index iq designates the color state of (i-quark; 

the low index iq designates the color state of g-quark. 
The contributions of the Feynman diagrams into the total amplitude can be written as follows 
(cg and Cc are electric charges of the q- and c-quarks correspondingly; the color coefficients are 
put into the braces): 



Ti = e, ■ n(g,)f'(p„ fe, 6„ J,,)vg, ■ {zZ-^-^^'^r.^Cl-}' (1°) 

= e, ■ ug,t\pg, fe, e„ J.Mq,) ■ {^/"^"^^'^Ci-C. J' (11) 

T, = eg- u{q,)t'{pg, qgg, eg, Jgg)v-^ ■ {ir'^'^nZ^r.X}, (12) 

T4 = e, ■ u,,t\pg, qgg, eg, Jgg)v{l^) ■ {^/"^"^^'^C^SJ, (13) 

T5 = ee ■ u,grv{qMQ,ho.Vg,/kl,, ■ (14) 



8 



Tio 
Til 

Tl2 
Tl3 
Ti4 
Tl5 
Tl6 

Tit 

Tl8 
Tl9 
^20 
T21 
T22 
223 
724 



= e„ 



uilgh'^^qgHlchaVcj/k 

Uqgl'^V {qq)Uc-yJaV {qc) / k 
u{(lq)l°'VqgUc^^aV{q^)/k 

■ {^ill^lic^iqiil}^ 

^ IqL I Lq Ic'^c ' 
' {^ial^lia^icic}^ 



n +3g +n 

qlq^ 

I ''/j - 

2 _ rj.n j.n ijg \ 
g,qq ' X'^icic iql liqi ^ 

2 

g,qq ^ lcac ig'- ".9- 
n ^jg 

c^c liq-* 



,2 



^c,qq^g^cy 
'^cgJqq^cr) 



'^cg^^'^c^qq 



'^c,qq^"/^cg 
Uc-yJqqVcg 



^c-f^^g^cqq 



^q,cc^g'^q'y 
'^qg J cc^q^ 



^qg^^^q,cc 



'^q,cc^^'^qg 
^57 J cc^qg 



^q-t^g^q,cc 



(15) 
(16) 
(17) 

(18) 
(19) 
(20) 
(21) 
(22) 
(23) 
(24) 
(25) 
(26) 
(27) 
(28) 
(29) 
(30) 
(31) 
(32) 
(33) 



The spinor states with the fixed spin projection on the axis z have been chosen in the capacity of 
two independent spinor states. In our calculations, the Dirac representation for the 7-matrices 
has been used. The spinor states can be written down as follows: 



u{p, +) 



VE + m 



uip, -) 



I 



vip, +) 



VeT 



m 



v{p, -) = 



VE + m 



( E + m \ 


Pz 

\Px + iPy J 

( \ 

E + m 
Px - iPy 
-Pz 

f Pz \ 

Px + iPy 


V E + m ) 

( Px- Wy \ 

Pz 



\ E + m J 



(34) 
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The 7-matrices in the Dirac representation are given by 
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1 
-1 



7 



a 
-a 



The Gell-Mann matrices have the form 



1 
-1 
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1 
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0\ 




' 1 
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V 1 





0. 
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v ^ 
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Vo 
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0/ 



/ \ 

-i 
\ i / 



A8 = ^ 



1 



/ 1 \ 

1 
\ -2 I 



and the ^-matrices have been chosen as f = |A*. 

Let us remind of the antisymmetric constant j^hc values, which have been required for our 
calculations 

/l23 = 1; 

/l47 = — /l56 = /246 = /257 = /345 = " /sGT = \i 
/458 = /678 = 

B. The soft process of c- and q-quark fusion into the (cg)-quarkonium 

To describe the fusion of c- and g-quarks into the (cg)-quarkonium we suppose that there are 
terms in the partonic distribution of D*-meson, which correspond to the valence quarks. We 
also suppose that the c- and g-quarks produced in the hard process transform into the valence 
quarks of the meson. The valence quark distributions in the system of infinite momentum have 
the following form: 

fc(^^P±) = fc{x,P±) - fc{x,P±), 
fg(x,p±) = fq{x,p±) - fq{x,p±). 

The averaged momentum fraction carried out by the valence quarks are 



(x^) = j (fp±dxx-f^{x,p±) 
(4) = j d^P±dxx ■ fg{x,p±) 



A 

mr,* 



(35) 
(36) 
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where A is the quark binding energy inside the meson. 
For (a;^) and (xp the following equation is valid 

(x:) + (xp ^ 1. (37) 

In our calculation we neglect the dispersion of the momentum fractions carried out by the 
quarks, or in other words, we consider ( p5| , p6| , p7D as absolutely precise equations. Also we 



suppose that the light quark mass is about A . 

In the doubly heavy quarkonium the quark pair is in the singlet state, because the octet 
one is suppressed by the third power of the relative velocity of quarks In the case of 

(cg)-quarkonium this suppression does not exist, and one has to take into account the both 
color states. 

The operators and (0(8)) describing the hadronization of singlet and octet quark pairs 

into the meson are nonperturbative, because they include the interaction on the scale about 
Aqq£)- III thfi framework of nonrelativistic potential model these operators correspond to the 
wave functions squared at the origin: (0(i_8))Inr ~ 1^(0) They are defined by 

(0(1)) = (^-g"' + ^) X (D*(p)|(c7,g)(g7.c)|D*(p)), (38) 

(0(8)) = ^ [-g^'" + j^J X {D*{p)\{c^,X^q){q^A'c)\D*{p)) — . (39) 

The probability W{c —>■ D*) of c-quark fragmentation into the D*-meson in e^e^-annihilation 
can be expressed through the operators (0(i)) and (0(8)) of our model 







where (O'^^ifiR)) = (0(i)) + i(0(8)), and 



J(r) = — 
^ ^ 27 



24 + 109r - 126r2 - IJAr^ - 89r^ r(7 - 4r - 3r^ + lOr^ + 2r^) , 

h — ^ In r 

15(1 -r)5 (l-r)6 



(41) 
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Fig. 1. The leading order QCD diagrams for the cg-state production in the 
g'7*-interaction. 
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Fig. 2. The D*-mcson production distribution versus Pt calculated in the 
framework of model under consideration (upper curve) in compar- 
ison with the fragmentation model prediction (lower curve) for the 
gf7-interaction at 200 GeV. 
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Fig. 3. The differential production cross sections of the £)*-meson for the 
deep inelastic ep-scattering: over the photon virtuality GeV^(a), 
the Bjorken x (b), the invariant mass of the final hadrons (c), the 
transverse momentum (d), the pseudorapidity (e) and the Feynman 
variable (/) in comparison with the experimental data of ZEUS 
Collaboration. 
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Fig. 4 The differential production cross sections of the £)*-meson for the 
deep inelastic ep-scattering over the invariant mass of the final 
hadrons (a), the transverse momentum (6), the Bjorken x (c), the 
pseudorapidity (d), the photon virtuality GeV^(e) and z{D*) in 
comparison with the experimental data of HI Collaboration. 
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Fig. 5. The experimental values of F2 (the ZEUS Collaboration) in com- 
parison with the discussed model predictions. 
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Fig. 6. The experimental values of F2 (the HI Collaboration) in compari- 
son with the discussed model predictions. 
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Fig. 7. The QCD leading order diagrams for the cg-state 
production in the ^7*-interaction. 
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Fig. 8 The singlet c^-pair production in the ^7— subprocess for the kine- 
matic region of HI Collaboration in comparison with the octet cq- 
pair production in the gi7-subprocess (the distributions are over the 
same variables as in Fig. 4) . 
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Fig. 9 The sum of singlet c^-pair production in the ^7— subprocess and 
the singlet cg-pair production in the g'7-subprocess in comparison 
with the data of HI Collaboration (the distributions are over the 
same variables as in Fig. 4) . 
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